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Abstract. A search is presented for a high mass neutral particle that decays directly to the e /i T final 
state. The data sample was recorded by the ATLAS detector in yfs = 7 TeV pp collisions at the LHC from 
March to June 2011 and corresponds to an integrated luminosity of 1.07 fb -1 . The data are found to be 
consistent with the Standard Model background. The high e /i^ mass region is used to set 95% confidence 
level upper limits on the production of two possible new physics processes: tau sneutrinos in an _R-parity 
violating supersymmetric model and Z'-\ike vector bosons in a lepton flavor violating model. 



Short-lived particles that decay into two oppositely 
signed leptons of different flavors, e /i^ (e/i), e ± r =F (er), 
or n (/ir), are predicted by a number of extensions to 
the Standard Model (SM). Examples include sneutrinos in 
i?-parity violating (RPV) supersymmetric (SUSY) mod- 
els jT] , and extra gauge Z' bosons with lepton flavor violat- 
ing (LFV) interactions [2] . This Letter reports a search for 
an excess of high invariant mass e/i (m ep ) events over SM 
predictions in pp collisions at y/s = 7TeV at the LHC. The 
e/i final state is chosen due to its clean detector signature 
and low SM background in the high m e)l region. Similar 
searches with the e/i final state have been reported previ- 
ously by the CDF, DO and ATLAS Collaborations 
[SHIS]. I n this Letter, we report an updated search with a 
data sample approximately 30 times larger than used for 
the previous ATLAS search [8] with improved sensitivity 
to new physics. 

The ATLAS detector [5] is a multi-purpose particle 
physics apparatus with a forward-backward symmetric cylin- 
drical geometry and near 4tt coverage in solid angl^E The 
inner tracking detector (ID) consists of a silicon pixel 
detector, a silicon microstrip detector, and a transition 
radiation tracker. The ID is surrounded by a thin su- 
perconducting solenoid providing a 2 T magnetic field 
and by a hermetic calorimeter system, which provides 
three-dimensional reconstruction of particle showers up to 

< 4.9. For \rj\ < 2.5, the electromagnetic calorime- 
ter is finely segmented and plays an important role in 



1 ATLAS uses a right-handed coordinate system with its ori- 
gin at the nominal interaction point (IP) in the centre of the 
detector and the z-axis along the beam pipe. The x-axis points 
from the IP to the centre of the LHC ring, and the y axis points 
upward. Cylindrical coordinates (R, <j>) are used in the trans- 
verse plane, <j> being the azimuthal angle around the beam pipe. 
The pseudorapidity is defined in terms of the polar angle 6 as 
r\ = -lntan(0/2). 



electron identification. The muon spectrometer (MS) is 
based on three large superconducting toroids arranged 
with an eight-fold azimuthal coil symmetry around the 
calorimeters. Three stations of drift tubes and cathode 
strip chambers enable precise muon track measurements, 
and resistive-plate and thin-gap chambers provide muon 
triggering capability. 

The data sample used in this analysis was collected us- 
ing single lepton (e or /i) triggers, between March and June 
2011. The total integrated luminosity is 1.07±0.04 fb -1 [TO]. 
The trigger efficiency is measured to be 100%, with a pre- 
cision of 1%, for e/i candidates that pass the default se- 
lection criteria described below. 

To select e/i candidates, the electron candidate is re- 
quired to have pr > 25GeV and to have pseudorapidity 
\ri\ < 1.37 or 1.52 < |?7| < 2.47. It is further required to 
pass the "medium" [11] quality definition, which is based 
on the calorimeter shower shape, track quality, and track 
matching with the calorimeter cluster. In addition, the 
electron is required to be isolated in the calorimeter with 
e AR<oa < 10GeVj where e£ r<oa is defined as the trans- 
verse energy deposited in the calorimeter within a cone of 
radius AR = y/ Arj 2 + A<fi 2 = 0.4 around the electron 
cluster. Corrections have been applied to account for en- 
ergy leakage from the electron and energy deposition in- 
side the isolation cone due to additional pp collisions. The 
muon candidate must be reconstructed in both the ID and 
the MS, and have p T > 25GeV and \rj\ < 2.4. Further- 
more, the muon is required to be isolated in the ID with 
p AR<0A < 10GeV; where p AR<0A ig defined as the scalar 

sum of the px of tracks associated to the primary vertex, 
within a cone of radius AR = 0.4 around the muon track. 
Only tracks with px > lGeV are used. Furthermore, only 
electrons separated from muons by AR > 0.2 are consid- 
ered. 
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The e/j, candidate events are required to have exactly 
one electron and one muon with opposite charge satisfying 
the above selection criteria. Furthermore, events have to 
contain at least one primary vertex reconstructed with at 
least three associated tracks of pr > 500 MeV. 

The SM processes that can produce an e/i signature 
can be divided into two categories: processes such as Z/-f* 
— > tt, tt, single top, WW, WZ and ZZ, which can pro- 
duce electrons and muons in the final state, and processes, 
referred to as fake background in this Letter, such as 
W/Z+'y, W/Z+jets and multijct events where the pho- 
ton or one or two jets are reconstructed as leptons. 

The contributions from processes listed in the first 
category as well as photon-related backgrounds are es- 
timated using Monte Carlo (MC) samples generated at 
y/s = 7TeV. The detector response simulation [T2] is based 
on the GEANT4 program [13j . Lepton reconstruction and 
identification efficiencies, energy scales and resolutions in 
the MC are corrected to the corresponding values mea- 
sured in the data in order to improve the modeling of the 
background. The MC predictions are normalized to the 
data sample based on the integrated luminosity and cross 
sections of various physics processes. Top production is 
generated with MC@NLO [T3] for ti and single top, the 
Drell-Yan process is generated with pythia [15], and the 
diboson processes are generated with herwig [TB] . Higher 
order corrections have been applied to the cross sections 
predicted by these generators [T71H81H9] . The W/Z + 7 
contribution in the fake background comes from the W{— > 
[iv)^i and Z{— > /x/z)7 processes, where the photon is re- 
constructed as an electron. This background is estimated 
using events generated with MADGRAPH )20j . 

The uncertainties for the ti and single top cross sec- 
tions are taken to be 10% [2T] and 9% [25], respectively. 
The cross sections for W/Z + 7, Z/7* ^ tt, WW, WZ 
and ZZ are assigned uncertainties of 10%, 5%, 7%, 7%, 
and 5%, respectively; these uncertainties arise from the 
choice of PDF, from factorisation and renormalisation scale 
dependence and from a s variations. The integrated lumi- 
nosity uncertainty and other smaller systematic uncertain- 
ties from the lepton trigger, reconstruction and identifica- 
tion efficiencies, energy (momentum) scale and resolution 
have been added in quadrature and are included in the 
total uncertainty. 

The remaining fake backgrounds arise from the W/Z+\ets cal1 ^ sl S mnc ant 
and multijet processes, where leptons are present from 
b- or c-hadron decays or at least one jet is misidentified 
as a lepton. Such lepton candidates are collectively re- 
ferred to as "non-prompt leptons" in this Letter. These 
jet fake backgrounds account for ^30% of the expected 
e/i data yield and are estimated from data using a 4 x 
4 matrix background estimation method described be- 
low. A looser lepton quality selection (called loose lepton 
here) is defined for each lepton type in addition to the 
default quality selection (called tight lepton here). For 
loose muons, the isolation requirement is dropped. For 
loose electrons, the "loose" electron identification criteria 
as defined in Ref. \TT\ are used and the isolation require- 
ment is also dropped. The tight and loose lepton selec- 



tions are then used to classify events where both leptons 
pass the loose requirements into four categories, depend- 
ing on whether both leptons subsequently pass the tight 
requirement (N pp ), only one lepton fails the tight require- 
ment and the other lepton passes the tight requirement 
(N p f or Nf p ), or both leptons fail the tight requirement 
(Nff). The sample composition can be estimated by solv- 



ing a linear system of equations: (N pp , N p f, Nf p , 



.N e f p f) T , where N efl (or N e 



Nff) T = 
,t) is the 



number of events with two prompt leptons (or two non- 
prompt leptons), while iV e „t and iV e t M are the numbers of 
events with one prompt lepton and one non-prompt lep- 
ton. The matrix e contains the probabilities for a loose 
quality lepton to pass the tight quality selection for both 
prompt and non-prompt leptons. The probability for prompt 
leptons (non-prompt leptons) is estimated by applying the 
loose and tight selections on Z/j* — > ee/[i[i events (a sam- 
ple of dijet events). To take into account the lepton px 
dependence of the two probabilities, the matrix equation 
is inverted for each event, giving four weights, correspond- 
ing to the four combinations of prompt and non-prompt 
leptons. These weights are then summed over all events 
to yield the total number of events with one or more non- 
prompt leptons. The overall jet fake background is found 
to be 1175 ± 32 (stat) events. The breakdown of these 
contributions is estimated to be N ep -t = 375 ± 30 (stat), 
N eip = 89 ± 13 (stat) and AT et|t t = 711 ± 8 (stat). The 
overall systematic uncertainty of 10% comes mainly from 
the uncertainty on the probability for a loose quality non- 
prompt muon to pass the tight quality selection. 

Table [1] shows the number of events selected in data 
and the estimated background contributions with their 
uncertainties (both statistical and systematic uncertain- 
ties are included). A total of 4053 e/i candidates are ob- 
served, while the expectation from SM processes is 4145 ± 
250 events. The m etl distribution is presented in Fig.[T]for 
data and background contributions. The distribution of 
observed events is compared to the expected background 
using a Kolmogorov-Smirnov test with statistical uncer- 
tainties only [23]. The test probability is 56%, consistent 
with the absence of a new physics signal. 

Table [2J shows the numbers of observed and predicted 
background events in eleven high e/i mass regions. Good 
agreement is found for all mass regions and no statisti- 
data excess is observed. Limits are set 
on the contributions of new physics processes to the high 
mass region from two scenarios: the production of v T in 
an RPV SUSY model and of an LFV Z 1 in extra-gauge 
boson models. 

The process dd — > v r — > e/i in a SUSY RPV model is 
considered. The RPV sneutrino couplings allowed in the 
supersymmetric Lagrangian are ,1. ,1 1. ■ \ t _, !.,(>, I )i . 
where L and Q are the lepton and quark SU (2) doublet su- 
perfields, and E and D denote the singlet fields for charged 
leptons and down type quarks, respectively. The indices 
i,j, k = 1,2,3 refer to the fermion generation numbers. 
The coupling constants A satisfy Xijk — ~^jik- Only the 
tau sneutrino is considered in this Letter since stringent 
limits already exist on the electron sneutrino and muon 
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Table 1. Estimated backgrounds in the selected sample, to- 
gether with the observed event yield. The total integrated lu- 
minosity is 1.07 fb _1 . 



Process 


Number of events 


tt 


1580 ± 170 


Jet fake 


1175 ± 120 


Z/'y* — > tt 


750 ± 60 


WW 


380 ± 31 


Single top 


154 ± 16 


W/Z + J 


82 ± 13 


wz 


22.4 ± 2.3 


zz 


2.48 ± 0.26 


Total background 


4145 ± 250 


Data 


4053 



Table 2. Estimated total backgrounds in the selected sample, 
together with the observed event yields for 11 high e/i mass 
regions. 





Data 


SM prediction 


> 200 GeV 


286 


288 ± 22 


> 250 GeV 


152 


136 ± 11 


> 300 GeV 


70 


67 ± 6 


> 350 GeV 


35 


34.0 ± 3.0 


> 400 GeV 


22 


17.7 ± 1.7 


> 450 GeV 


10 


10.5 ± 1.2 


> 500 GeV 


7 


6.8 ± 0.9 


> 550 GeV 


3 


4.3 ± 0.6 


> 600 GeV 


3 


2.4 ± 0.4 


> 650 GeV 


1 


1.49 ± 0.31 


> 700 GeV 





1.07 ± 0.25 



sneutrino [TJ. By fixing all RPV couplings except A 311 (i> r 
to dd) and A312 (p T to e/i) to zero, and assuming that 
v T is the lightest supersymmetric particle, the contribu- 
tions to the e/i final state originate from the v T only. The 
cross section is 0.154 pb for mp T = 650GeV, A 311 = 0.10 
and A312 = A321 = 0.05 [24]. The total decay width is 
r„ T — (3A 3 2 n + 2A3 12 )toj> t /167t. Using couplings that are 
consistent with the current limits, the decay width is less 
than lGeV for m,v T — ITeV, which is well below the con- 
tribution from detector resolution. MC samples with v T 
masses ranging from 0.1 to 2TeV are generated with her- 
wig [TBpB] . 

An e/i resonance also appears in models containing 
a heavy neutral gauge boson, Z' [26] . with non-diagonal 
lepton flavor couplings. Rare muon decay searches have 
placed extremely stringent limits on the combination of 
the mass and the coupling to ee and e/i in such models [5] . 
The e/i searches at hadron colliders are not able to match 
the sensitivity of dedicated /i —¥ e conversion experiments. 
A limit on the production cross section times branching ra- 
tio to e/i is placed on the Z'-like boson model to represent 
the production of vector particles that can decay to the 
e/i final state. To calculate the efficiency and acceptance, 
the Z' is assumed to have the same quark and lepton cou- 
plings as the SM Z except a non-zero Z' to e/i coupling, 



which is assumed to be the same as the Z' to ee coupling. 
The cross section is 0.61 pb for m z > = 700GeV [27]. MC 
samples with Z' masses ranging from 0.7 to 2 TeV are 
generated with pythia. 

Both v T and Z' samples are processed through the 
standard chain of the ATLAS simulation and reconstruc- 
tion. The overall product of acceptance and efficiency is 
36% for mj> 7 = lOOGeV and increases to 64% for = 
ITeV. The corresponding number is ~ 60% for Z' with 
mass mz' = 700GeV to mz' = 2TeV. The predicted m efJi 
distributions for a v T with mp r = 650GeV and a Z' with 
mz> = 700GeV are also shown in Fig. [TJ 



> 
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Fig. 1. Observed and predicted e/t invariant mass distribu- 
tions. Signal simulations are shown for m^ T — 650GeV and 
m z , = 700GeV. The couplings X' 3U = 0.10 and A312 = 0.05 
are used for the RPV v T model. The production cross section 
is assumed to be the current published limit of 0.178 pb for 
the LFV Z' model [8]. The ratio plot at the bottom includes 
only statistical uncertainties. 

The rn eAI spectrum is examined for the presence of a 
new heavy particle. For each assumed ma r value in the 
range 100 GeV to 2 TeV, a search region, which depends 
on the simulated e/i mass resolution, is used3- The number 
of observed and predicted background and signal events 
in each search range are used to set an upper limit on 
<j(pp —> v T ) x BR(£ r — > e/i). A Bayesian method [25] is 
used with a uniform prior for the signal cross section for a 
given TO£ x . Fig. shows the expected and observed 95% 
confidence level (C.L.) limits, as a function of m„ T , to- 
gether with the limits previously published by ATLAS [8J , 
which were based on 35 pb^ 1 of data, and the expected 
±1 and ±2 standard deviation uncertainty bands. For a 
v T with a mass of 100 GeV (1 TeV), the limit on the 
cross section times branching ratio is 135 (4.5) fb. The 
limits obtained extend 7 (34) times beyond the previous 

2 The search region is normally defined to be (m^ — 3a, 
ma T + 3a), where a is the expected m efl resolution (e.g., a — 
11 GeV for = 400GeV). If mp T - 3a < 700GeV and 

ma T +3<7 > 700GeV, the region above mo T —3a is used. If — 
3a > 700GeV, the region above 700GeV is used. The mass 
window changes around 700GeV because the MC statistics is 
not sufficient in the m efl > 700GeV region. 
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ATLAS results. The theoretical cross sections for A 311 = 
0.10, A 3 i 2 = 0.05 and X' 3n = 0.11, A 3 i 2 = 0.07 are also 
shown. Tau sneutrinos with a mass below 1.32 (1.45)TeV 
are excluded, assuming coupling values A 311 = 0.10 and 
A312 = 0.05 (X' 3n = 0.11 and A 3 i 2 = 0.07). The limits 
are significantly better than the limits from the previous 
ATLAS analysis using 35 pb" 1 of data. The 95% C.L. 
observed upper limits on A 311 as a function of mp T are 
shown in Fig. [2}d for three values of A 3 i2, together with 
the exclusion region obtained from the DO experiment [7] 
and previously by the ATLAS experiment [5]. The limits 
on A 311 are tighter than the DO results for m,v r > 270GeV 
sneutrinos assuming A 31 2 = 0.07. Better sensitivity can be 
obtained for m^ x < 270 GeV by applying selection cuts 
on missing transverse energy and number of jets in the 
event to improve the signal and background ratio, but it 
will make the search model-dependent. 

A similar method is used to set limits on the LFV Z'- 
like vector boson; however, as opposed to the sneutrino 
limits, a unique mass window is defined for each poten- 
tial signal mass. The 95% C.L. upper limits on a(pp — > 
Z') x BR(Z' — > e/i) are shown in Fig. [3] The expected 
limit is the same as the observed limit for the high mass 
points because both the median background event count 
expectation and the observed number of events are zero. 
For a Z' with mass of 0.7 TeV (1.0 TeV), the limit on the 
cross section times branching ratio is 9.6 fb (4.8 fb). This 
result improves upon previous ATLAS limits by roughly 
a factor of 20 (40). 

In conclusion, a search has been performed for high 
mass efj, events using pp collision data at s/s = 7TeV 
recorded by the ATLAS detector. The observed m ejl distri- 
bution is found to be consistent with SM predictions. With 
no evidence for new physics, 95% C.L. exclusion limits 
are placed on the production cross sections and RPV cou- 
pling values of the tau sneutrinos in an RPV SUSY model, 
and tau sneutrinos with a mass below 1.32 (1.45)TeV 
are excluded, assuming coupling values A 311 = 0.10 and 
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A 3 i2 = 0.05 (A 311 = 0.11 and A. 



312 = 0.07). The results 
presented here are the most stringent results to date for 
m,y T > 270 GeV. More stringent constraints are also set 
on the production cross sections of Z' bosons in an LFV 
model. These two benchmark models can be used to rep- 
resent the production of any narrow scalar and vector par- 
ticles that can decay to the e/j, final state. 
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A312 = 0.05 and A 3 n = 0.11, A312 = 0.07 calculated using 
MADGRAPH with next-to-leading order fc-factors applied [21] 
are also shown, (b) The 95% C.L. upper limits on the A311 cou- 
pling as a function of mp T for three values of A 3 i2- The regions 
above the three curves represent ranges of A311 values that 
are excluded. These results are compared with the exclusion 
regions obtained from the DO experiment and the previously 
published ATLAS analysis. The cross section times branching 
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which causes the weak dependence of the A 311 limits on A312 
for low mass tau sneutrinos. 
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Fig. 3. The observed 95% C.L. upper limits on a(pp — > Z') x 
BR(Z' — » efi). The expected limits are also shown together 
with the expected ±1 and ±2 standard deviation uncertainty 
bands. The observed and expected limits overlap as discussed 
in the text. 
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